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Introduction 

The  long-term  objective  of  this  effort  is  to  develop  a  generic  gene  shuffling-based  technology 
to  rapidly  screen  libraries  of  1010  proteins/peptides  encoded  by  DNA  libraries,  for  identifying 
biomolecules  that  can  intercept  both  existing  and  emerging  organophosphate-based  chemical 
warfare  nerve  agents  (CWNAs).  Enzymes  identified  in  these  screens  should  be  capable  of 
catalytically  neutralizing  the  target  agent  under  physiological  conditions,  thereby  providing  a 
basis  for  development  of  a  new  generation  of  therapeutic  agents  against  CWNAs.  The  major 
milestone  was  to  integrate  established  components  of  enhanced  molecular  evolution  techniques 
so  as  to  provide  a  means  of  miniaturizing  existing  low-throughput  assays,  thereby  dramatically 
increasing  both  sensitivity  and  throughput.  The  expressed  biomolecules  were  allowed  to 
intercept  the  CWNA,  preventing  interaction  with  its  target  (viz,  AChE),  and  genes  code  for 
effective  interceptors  were  isolated.  Uniquely,  this  screen  is  for  directly  detoxifying  the  CWNA, 
not  simply  for  binding  it,  thus  allowing  identification  of  biomolecules  that  prevent  its  action  by 
degradation  at  the  desired  rate.  This  technology  was  envisaged  to  provide  rapid  discovery  of 
pretreatment  and  post-challenge  therapeutic  drugs  against  existing  and  emerging  CWNA  threats, 
and  thus  to  shorten  the  time  from  emergence  of  a  threat  to  identification  of  potential  counter¬ 
measures  to  a  few  days  or  weeks. 

G-  and  V-type  type  nerve  agents:  tabun  (GA),  sarin  (GB),  soman  (GD),  cyclosarin  (GF)  and 
VX  (Figure  la),  are  highly  toxic  organophosphates  (OPs)  that  penetrate  the  body  both  by 
inhalation  and  skin  absorption4.  They  contain  a  chiral  phosphorous  center,  whose  optical  isomers 
differ  in  their  toxicity.  The  Sp  isomers  of  GB,  GD,  and  GF  inhibit  acetylcholinesterase  (AChE) 
>1000  fold  faster  than  their  respective  Rp  isomers  and  are  considered  to  be  much  more  toxic5. 
The  Sp  isomer  of  VX  inhibits  AChE  approximately  100-fold  faster  than  the  Rp  enantiomer.  For 
GD,  in  addition  to  its  chiral  phosphorous,  a  chiral  carbon  center  (Figure  la)  produces  4 
diastereo-isomers:  SpSc,  SpRc,  RpSc  and  RPRC.  However,  since  the  chirality  of  the  phosphorous  in 
GD  has  a  much  greater  influence  on  AChE  toxicity  than  the  chirality  of  its  alpha  carbon6,  only 
two  of  its  isomers  are  considered  toxic:  SpSc  and  SpRc-  In  the  case  of  GA,  which  contains  a 
cyanide  in  place  of  a  fluoride  (Figure  la),  nomenclature  conventions  dictate  an  Rp  configuration 
for  its  toxic  isomer.  Following  their  entrance,  G-agents  bind  to  and  inhibit  AChE,  thus  causing  a 
cholinergic  crisis  that  may  result  in  respiratory  failure  and  death7.  Current  treatments  of  OP 
intoxication  are  symptomatic  and  do  not  act  to  reduce  blood  OP  levels  and,  thereby,  to  prevent 
AChE  inhibition  ’  .  A  better  treatment  strategy  might  be  achieved  by  employing  OP  hydrolyzing 
enzymes  as  prophylactic  drugs  that  would  promote  rapid  OP  degradation  in-vivo10'12. 
Unfortunately,  most  OP-hydrolyzing  enzymes  isolated  thus  far  show  preference  towards  the  non¬ 
toxic  isomers  of  nerve  agents,  and  hydrolyze  the  toxic  isomers  with  low  catalytic 
efficiencies13’14.  Until  this  project  was  initiated,  there  were  no  available  enzymes  capable  of 
hydrolyzing  the  toxic  isomers  of  nerve  agents  at  a  rate  that  would  permit  them  to  offer  effective 
prophylactic  protection  at  low  protein  doses  (<  50mg/70kg).  Notably,  treatment  with  the  most 
advanced  non-catalytic  enzyme  scavenger,  human  butyrylcholinesterase,  requires  the 
administration  of  a  dose  of  several  hundred  milligrams  per  person  in  order  to  confer  reasonable 
protection  against  nerve  agents15. 

PON1  is  a  mammalian  lactonase  that  exhibits  a  promiscuous  OP-hydrolyzing  activity3'16.  It 
resides  on  cholesterol-carrying  HDL  blood  particles  where  it  plays  a  role  in  drug  metabolism  and 
in  the  prevention  of  atherosclerosis17’19.  Human  PON1,  and  the  other  mammalian  PONs, 
hydrolyze  G-agents  at  low  catalytic  efficiencies  and  react  primarily  with  their  less  toxic  Rp 
isomers'  .  Using  directed  evolution  we  obtained  a  recombinant  PON1  (rePONl)  variant,  G3C9, 
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which  expressed  well  in  soluble  form  in  E  coli,  that  was  selected  as  the  wt  lead  variant  for 
directed  evolution. 

Over  the  last  5  years  of  performance  in  the  framework  of  this  project,  we  have  demonstrated 
the  potential  of  directed  evolution,  combining  random  and  designed  mutations  based  on  3D 
structures,  to  generate  mutants  of  a  recombinant  mammalian  PON1  (rePONl)  with  catalytic 
proficiency  well  above  that  of  the  wild  type  '  .  We  have  evolved  mutants  that  hydrolyze  the  toxic 
isomers  of  the  coumarin-based  nerve  agent  surrogates  >105-fold  faster  than  rePONl  and  human 
PON1.  We  have  also  developed  a  safe  in  situ  protocol  for  generating  the  corresponding 
fluoridates  ( i.e the  threat  G  agents),  and  showed  that  kcat/Kra  values  for  these  variants  are  as 
high  as  2-5xl07  M ’min"1,  values  approaching  the  theoretical  estimated  minimal  requirement  for 
efficient  prophylactic  protection  at  reasonable  doses  (~50  mg/70  kg)  without  the  need  for  post¬ 
exposure  therapy.  In  addition,  we  have  shown  that  PON1  variants  evolved  by  selection  with  VX 
are  capable  of  hydrolyzing  V-type  agents,  and  anticipate  that  the  current  libraries  will  be  further 
developed  by  use  of  directed  evolution  to  evolve  one  or  more  variants  that  will  be  qualifie  as  V- 
type  agent  antidotes.  It  seems  likely  that  not  more  than  two  PON1  variants  will  be  required  to 
provide  broad-spectrum  protection  against  all  currently  prevalent  nerve  agents. 


Body 

I.  Specific  Aims 

1.  The  development  of  high-throughput-assays  for  OP  hydrolase  variants  exhibiting  high 
specificity  factors  and  turnover. 

2.  Provision  of  proof-of-concept  for  the  proposed  core  technology  employing  directed 
evolution  of  new  recombinant  PON  and  AChE  variants. 

3.  Isolation  of  interceptors  for  G-  and  V-type  nerve  agents,  and  expression  in  soluble  form. 

4.  Design,  generation  and  selection  of  2nd  generation  libraries  for  V-  and  G-type  agents. 

5.  Large-scale  production  of  selected  enzyme  candidates,  and  their  kinetic,  structural  and 
pharmacological  evaluation 

6.  Establishment  of  “off-the-shelf’  libraries  for  rapid  identification  of  antidotes  against 
emerging  future  threats 


II.  Significance  to  the  goals  of  Counter  ACT 

The  proposed  approach  opens  new  opportunities  for  rapid  identification,  characterization  and 
implementation  of  novel  countermeasures  against  CW  agents.  It  will  significantly  decrease  the 
time  interval  between  the  appearance  of  a  new  threat  and  the  discovery  of  potential  antidotes  to 
counteract  it.  The  major  benefits  will  be  one  or  more  products  capable  of  efficient  catalytic 
hydrolysis  of  G-  and  V-type  nerve  agents,  as  well  as  of  gene  libraries  derived  from  existing 
enzymes  that  can  be  used  “off-the-shelf’  to  isolate  new  protein  variants  for  almost  any  nerve 
agent  or  toxic  industrial  chemical  serving  as  a  target  for  the  screen. 
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III.  Detailed  Tasks  -  Performance  and  accomplishments. 

A.  Generation  of  the  round  zero  library  of  G-type  CWNA  hydrolyzing  variants1: 

We  used  a  recombinant  variant  dubbed  rePONl-G3C9  (referred  to  as  rePONl  in  the 
following)  that,  unlike  human  PON1,  is  expressed  in  soluble  and  functional  form  in  E.  coli.  Its 
amino  acid  sequence  is  94%  and  85%  identical  to  rabbit  and  human  PONls,  respectively".  The 
enzymatic  specificity  of  rePON  1  is  essentially  identical  to  that  of  human  PON  1 ,  yet  its  stability 
is  much  improved,  and  its  in  vivo  applicability  and  efficacy  have  been  confirmed"  .  For  safety 
and  for  facile  detection,  we  initially  used  fluorogenic  coumarin  analogues  (Figure  lb),  including 
CMP-coumarin  that  inhibits  AChE  with  a  k,  that  is  20-40  fold  lower  than  cyclosarin  itself 
(depending  on  the  source  of  AChE)  . 

We  earlier  isolated  several  variants  of  rePONl  with  enhanced  activity  towards  a  racemic 
mixture  of  CMP-Coumarin  by  screening  ‘neutral  drift’  libraries  of  rePONl  (e.g.  1G3,  2G9)  . 
The  most  active  variant  was  3B3,  which  had  ~250-fold  higher  catalytic  efficiency  (kcat/Km 
20xl06  Mimin'1)  compared  to  the  wild-type-like  rePONl  (kcat/Km  0.08xl06  Mimin'1).  Although 
hydrolysis  by  rePONl-3B3  was  restricted  to  the  R p  isomer  (not  shown),  we  could  use  its  high 
catalytic  efficiency  and  Rp-  stereoselectivity  to  isolate  the  Sp  isomers  of  CMP-coumarin  and  IMP- 
coumarin  from  the  corresponding  racemates,  and  apply  them  in  subsequent  screens.  We  also 
identified  low  activity  of  H115W  and  V346A  rePONl  mutants  towards  Sp-CMP-coumarin20. 
However,  their  activity  with  Sp-CMP-coumari  n  was  too  low  for  detection  under  library  screening 
conditions.  We  therefore  used  for  the  initial  rounds  of  screening  the  coumarin  analog  of  sarin  ( 
IMP-coumarin),  a  less  bulky  G-agent  surrogate  whose  Sp  isomer  is  more  reactive  with  PON1 
(Table  1). 

Random  mutagenesis  of  rePON  1-H115W-V346A,  followed  by  screening  of  the  resulting 
library  in  96-well  plates  with  Sp-IMP-coumarin,  yielded  several  improved  variants  that  typically 
carried  one  mutation  in  addition  to  H115W  and  V346A.  A  second  round  of  mutagenesis  and 
screening  with  Sp-IMP-coumarin  led  to  isolation  of  variants  in  which  the  V346A  mutation  was 
removed  and  the  H115W  and  F222S  mutations  dominated.  As  the  evolving  variants  became 
more  reactive  towards  the  Sp  isomer,  the  3rd  generation  library  could  be  screened  with  both  Sp- 
IMP-  and  Sp-CMP-coumarin.  Indeed,  this  round  resulted  in  several  variants  with  improved 
activities  towards  Sp-CMP-coumarin  (e.g.,  3A7,  8C8;  Table  2).  However,  since  the  4th  round  of 
mutagenesis  and  screening  yielded  no  further  improvements,  we  designed  a  structure -based 
targeted  library  and  subjected  it  to  high-throughput  screening  (>106  variants  per  run)  by  FACS 
sorting23  24,  as  described  below. 

The  targeted  substitutions  library  was  based  on  PONl’s  active-site  structure.  In  particular,  a 
recently  obtained  crystal  structure  of  the  re-G3C9-Hl  15W  indicated  movements  of  several  side- 
chains  in  response  to  this  mutation,  including  those  of  residues  69  and  134,  that  are  in  direct 
contact  with  W115,  and  of  the  more  remote  residues  346,  347  and  348  (Ben-David  el  al. 
unpublished  data).  We  therefore  generated  a  library  by  randomizing  these  positions  and  those  of 
residues  115  and  222  that  had  been  found  to  be  mutated  in  all  active  variants  of  the  3ld  round 
(Table  3).  We  applied  an  oligo  spiking  strategy  that  incorporated  the  randomizing  oligos  onto  re¬ 
PON  1-H 11 5W  in  a  combinatorial  manner,  so  that  each  library  variant  carried  on  average  4 
mutated  positions25.  Due  to  the  intense  level  of  mutagenesis,  and  the  targeting  of  the  active  site, 
the  libraries  consisted  mostly  of  inactive  variants.  To  purge  inactive  library  clones,  we  employed 


1  For  experimental  details  and  additional  documented  results  please  refer  to:  Gupta  et  al  (201 1)  Nat  Chem  Biol,  7, 

120-125  (+  Supplementary  Material) 


-6- 


a  high-throughput  FACS  screen  using  a  fluorogenic  phosphotriester  dubbed  DEPCyC,  whose 
susceptibility  was  found  to  correlate  well  with  the  activity  on  Sp-CMP-coumarin"  .  E.  coli  cells 
transformed  with  the  plasmid  library  were  compartmentalized  in  water- in-oil  emulsion  droplets. 
The  fluorogenic  substrate  was  added,  and  the  primary  emulsion  droplets  were  converted  to 
double-emulsion  droplets  that  were  sorted  by  FACS.  Cells  in  isolated  droplets  were  plated, 
picked,  and  assayed  in  96-well  plates  for  Sp-CMP-coumarin  activity. 

The  most  active  variants  isolated  from  this  5th  round  carried  mutations  L69G/A  and  H134R, 
in  addition  to  H115W  and  F222S  that  had  appeared  in  earlier  rounds.  These  variants  were 
shuffled,  together  with  random  mutagenesis  at  low  rates  (~1.7  amino  acid  exchanges  per  gene), 
and  the  resulting  6th  round  library  was  sorted  by  FACS,  and  then  screened  in  96-well  plates.  The 
most  improved  variants  carried  five  key  mutations:  L69G,  H115W,  H134R,  F222S  and  T332S 
(Table  4).  The  best  variant,  3D8  exhibited  a  kcat/Km  value  of  1.2xl07  Mimin'1  with  Sp-CMP- 
coumarin.  Furthermore,  3D8  and  other  6th  round  variants  exhibited  similar  rates  with  both  the  Sp 
and  Rp  isomers  as  indicated  by  the  complete  hydrolysis  of  the  racemic  CMP-coumarin  with 
monophasic  kinetics  (Fig.  2). 

Along  with  selection  of  variants  displaying  higher  rates,  we  decreased  the  concentration  of 
the  CMP-coumarin  substrate  so  as  to  permit  isolation  of  variants  with  improved  Km  as  well  as 
kcat.  However,  to  be  able  to  screen  at  concentrations  that  correspond  to  the  very  low  toxic 
concentrations  of  cyclosarin  in  vivo  (-1  pM)415,  and  for  variants  that  efficiently  degrade 
cyclosarin  itself,  we  developed  a  screen  based  on  monitoring  the  rescue  of  AChE,  the  OP’s 
physiological  target.  AChE  was  added  to  crude  bacterial  lysates  expressing  the  library  PON1 
variants.  The  OP  was  added,  and  the  residual  activity  of  AChE  was  subsequently  measured  using 
a  chromogenic  assay  to  indicate  the  level  of  OP  degradation  by  the  tested  variant.  The  FACS 
sorted  6th  round  library  was  re-screened  in  96-well  plates  using  the  AChE  assay,  and  1  pM 
racemic  CMP-coumarin.  We  screened  730  randomly-picked  colonies,  and  isolated  13  variants 
that  were  improved  by  2-12  fold  relative  to  3D8. 

Although  at  this  stage  we  had  identified  variants  that  exhibited  sufficiently  high  catalytic 
proficiency,  these  had  been  selected  and  tested  with  coumarin  surrogates.  The  fluoride  leaving 
group  of  the  actual  threat  agents  differs  substantially  from  coumarin  (Fig.  1)  -  inasmuch  as 
fluoride  is  both  more  reactive  and  much  smaller.  However,  the  toxicity  of  nerve  agents  prevents 
their  use  in  ordinary  labs.  We  therefore  developed  a  non-hazardous  screening  protocol  based  on 
CMP-F  generated  in  situ,  in  dilute  aqueous  solution,  by  replacing  the  coumarin  leaving  group  of 
CMP-coumarin  with  fluoride.  This  exchange  was  spectroscopically  monitored  by  following  the 
release  of  the  coumarin  in  the  presence  of  NaF.  The  inhibition  of  To  AChE  by  CMP-F  thus 
generated  in  situ  proceeded  40-fold  faster  than  with  Sp-CMP-coumarin,  with  the  expected  k,  of 
1.3xl09  Mimin'1.  Using  the  in  situ  generated  CMP-F,  we  measured  the  kcat/Km  values  of  the 
evolved  variants  under  pseudo-first-order  conditions  (CMP-F  <50  nM,  well  under  a  likely  Km). 
Encouragingly,  the  activities  with  the  coumarin  and  fluoridate  Sp  isomers  were  comparable,  and 
at  least  three  variants  (0C9,  2D8,  1A4)  exhibited  kcat/Km  values  of  >107  Mimin'1  with  CMP-F 
(Table  5).  The  AChE  protection  assays  thus  confirmed  the  ability  of  the  evolved  variants  to 
protect  AChE  from  cyclosarin  in  vitro,  and  validated  the  coumarin  analogues  as  faithful 
surrogates  of  the  actual  G-agents.  The  assay  also  confirmed  that  the  starting  point,  rePONl- 
G3C9  is  much  more  active  towards  CMP-F  than  towards  CMP-coumarin  (Table  2),  as  is  human 
PON126,  although  the  kcat/Km  (-105  Mimin'1)  is  >  100-fold  too  low  for  in  vivo  detoxification 
using  reasonable  amounts  of  enzyme. 
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The  evolved  variants  were  sufficiently  active  to  permit  library  screens  using  the  in  situ 
generated  agent.  This  approach  is  highly  attractive,  since  the  assay  of  AChE  protection  against 
the  actual  threat  agent  mimics  the  in  vivo  protection  challenge  whereby  the  catalytic  scavenger 
must  be  sufficiently  active  to  intercept  the  threat  agent  before  the  latter  reacts  with  AChE.  We 
therefore  re-screened  the  13  most  improved  variants  from  the  last  round,  using  CMP-F  at  the 
expected  plasma  concentration  for  lxLDso  exposure  (1  pM).  Nine  variants  exhibited  improved 
activities  relative  to  3D8.  Of  these,  3  variants  (4E9,  5F3  and  6A3)  exhibited  the  highest  specific 
activity  upon  determination  of  the  amount  of  soluble  expressed  protein.  Following  sequencing 
and  protein  purification,  we  identified  4E9  as  the  most  active  variant  with  both  Sp-CMP- 
coumarin  and  Sp-CMP-F  (kcat/Km  =  2.23xl07  Mimin'1,  and  1.7xl07  Mimin'1,  respectively). 

B.  Generation  of  advanced  PONls  hydrolyzing  G-type  nerve  agents 

We  generated  gene  libraries  derived  from  the  lead  2D8  of  round  zero  above  (Table  4)  using 
structure-guided  targeted  mutagenesis.  The  3D  structure  of  2D8  copied  into  the  G2E6  scaffold 
(Fig.  3)  revealed  a  flexible  loop  that  covers  the  active  site.  Replacement  of  residues  that  were 
already  fixed  in  2D8  enabled  us  to  generate  smaller  and  more  focused  libraries  that  could  be 
screened  by  our  low-throughput  96-well  plate  screen.  Several  strategies  were  applied  for  the 
generation  of  advanced  mutant  libraries  beyond  round  zero,  as  detailed  below. 


Strategy  1:  Targeted  diversification  of  active -site  residues 

Examination  and  comparison  the  structures  of  wt  rePONl  (PDB  code  1V04  )  itself,  of 
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rePONl  in  complex  with  the  inhibitor  2-hydroxyquinone  (PDB  code  3SRG'  ),  and  of  a  docking 
model  of  rePONl  with  the  OP  pesticide  paraoxon  ,  pointed  to  active-site  positions  that  could 
affect  OP  binding  and  hydrolysis  (e.g.,  residues  70,  71,  196,  240  and  292).  We  also  assumed  that 
the  activity  of  evolved  variants  could  be  further  optimized  by  reexamining  the  mutagenesis  of 
active-site  residues  that  were  substituted  in  the  earlier  rounds  of  directed  evolution  towards  GF 
hydrolysis  (e.g.,  positions  69,  115,  134).  We  constructed  the  library  by  spiking  mutations  in  these 
8  active-site  positions  into  the  2D8  gene  in  a  combinatorial  manner  using  the  ISOR  method'  .  In 
general,  residues  were  mutated  to  amino  acids  with  similar  physico-chemical  properties  although 
more  drastic  changes  were  also  included  (e.g.  Gly69Leu/Val/Ile/Ser).  Sequencing  of  randomly 
selected  clones  revealed  that  the  unselected  library  contained  an  average  of  4  mutations  per 
clone,  with  each  clone  exhibiting  a  different  mutational  composition.  This  round  1  library  was 
then  screened  with  in-situ  generated  GD  and  GB,  using  the  AChE  inhibition  assay.  This  assay 
measures  the  ability  of  PON1  variants  to  prevent  loss  of  AChE  activity  by  rapidly  hydrolyzing 
the  toxic  isomer  of  the  added  OP. 

By  the  end  of  round  1  (see  Table  6a)  all  improved  variants  had  acquired  mutations  at  two 
positions,  69  and  115,  which  had  previously  been  mutated  to  improve  GF  hydrolysis1.  Upon 
selection  towards  a  broader  range  of  G-agents,  these  residues  changed  again.  In  addition, 
improved  variants  had  acquired  a  few  more  mutations,  primarily  at  positions  70,  71,  196,  240  or 
292.  The  changes  in  residues  69  and  115  led  us  to  explore  the  re-optimization  of  another  key 
active-site  residue,  222,  whose  mutation  to  Ser  had  led  to  increased  GF  hydrolysis1.  Residue 
Ser222  was  therefore  targeted  for  mutagenesis  in  the  2nd  round  library  that  was  screened  for  GB, 
GD  and  GF  neutralization  (Table  6b).  We  speculated  that  mutations  to  hydrophobic  residues 
(Feu,  Val,  He,  Met)  or  to  Cys  might  be  beneficial.  Indeed,  96%  (24/25)  of  improved  clones  from 
the  2nd  round  library  were  mutated  at  position  222,  mostly  to  Met  (17/24). 
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In  the  3rd  round  library  (Table  6b),  we  targeted  various  residues  whose  side-chains  are 
involved  in  structuring  the  top  of  PONl's  active-site  cleft.  These  included  residues  50,  197,  291, 
294,  346,  347  and  348.  In  addition,  we  retargeted  positions  196  and  292  that  had  been  explored 
during  the  generation  of  the  lst-round  library,  but  acquired  no  mutations.  On  average,  each  clone 
of  the  naive  (i.e.,  unselected)  3  round  library  contained  2+1  substitutions  at  the  targeted 
positions.  Following  screening  with  GB  and  GD,  we  obtained  18  improved  clones  in  which 
positions  197  and  291  were  primarily  mutated,  i.e.,  22%  (4/18)  and  44%  (8/18)  of  selected 
clones,  respectively.  Some  positions  did  not  accept  any  mutations  (196,  294,  346,  347,  348),  and 
others  were  mutated  in  single  clones  (50,  292). 


Strategy  2:  Random  shuffling  of  improved  variants 

Shuffling  of  selected  variants  can  combine  beneficial  mutations  that  were  acquired  in 
separate  variants,  or  eliminate  deleterious  mutations  from  individual  variants.  In  particular,  since 
the  mutational  diversity  in  our  libraries  was  large  compared  to  the  number  of  clones  screened 
(screening  covered  between  0.07  to  4xl0"4  of  the  theoretical  library  diversity),  most  active 
variants  isolated  exhibited  different  mutational  compositions.  In  fact,  shuffling  of  improved 
variants  was  applied  in  every  round  in  conjunction  with  the  ISOR  protocol  that  was  used  to 
introduce  the  targeted  mutations.  For  the  round  4  library,  however,  shuffling  was  applied  on  its 
own  with  no  spiking  of  targeted  mutations. 


Strategy  3:  Targeted  mutagenesis  of  the  flexible  active- site  loop 

PONl's  longest  active  site  loop  (residues  70-78)  is  highly  flexible,  and  is  not  observed  in  the 
crystal  structures  of  wild-type-like  rePONl  in  its  uncomplexed  form  .  However,  upon  binding  of 
the  lactam  inhibitor,  2-hydroxyquinoline  (2HQ),  the  loop  became  structured,  thus  completing  the 
active-site  cleft  in  the  crystal  structure”  .  Residues  Tyr71  and  Ile74,  in  particular,  comprise  part 
of  the  active-site  wall,  and  make  contact  with  the  2HQ.  In  addition,  the  loop  configuration,  and 
the  position  of  Tyr71  in  particular,  seem  to  vary  between  different  substrates  .  Indeed,  mutations 
in  positions  70  (Lys  to  Asn  or  Gin)  and  71  (Tyr  to  Phe  or  He)  appeared  in  many  improved  clones 
from  round  2.  During  round  3,  we  therefore  designed  a  library  devoted  to  examining  the  effects 
of  active-site  loop  mutations.  Specifically,  given  the  very  small  size  of  the  fluoride-leaving  group 
of  G-agents,  we  explored  loop  mutations  that  introduced  large  side-chains,  such  as  Trp,  Phe  and 
Tyr  that  might  reduce  the  active-site  volume  and  improve  substrate  binding.  We  also  attempted 
to  affect  the  loop's  configuration  by  mutating  its  residues  to  Gly  and  Pro.  Since  we  introduced 
amino-acid  substitutions  using  short  oligonucleotides  (21-29  bp),  most  substitutions  were 
mutually  exclusive,  and  each  clone  carried,  on  average,  only  one  loop  mutation. 


Strategy  4:  Introduction  of  ancestral  mutations. 

Ancestral  mutations,  i.e.,  substitutions  at  residues  that  have  appeared  along  the  evolutionary 
history  of  a  given  protein,  have  been  shown  to  be  powerful  modulators  of  a  protein's  stability  and 
function”  "  .  A  library  based  on  active-site  substitutions  borrowed  from  the  predicted  ancestors 
of  vertebrate  PONs  has  been  constructed  .  This  library  was  screened  for  OP-hydrolyzing 
activities  using  CMP-coumarin,  and  several  ancestral  mutations  that  were  able  to  increase  this 
activity  in  variants  were  identified”  .  We  included  these  mutations  (Leu55Ile,  Ile74Leu, 
Aspl36His,  Prol89Gly)  together  with  similar  ones  (55Met/Val,  136Gln,  189Ser)  in  the  round  3 
library,  at  an  average  frequency  of.  0.5  mutations  per  gene  in  the  unselected  library. 
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Table  6b  summarizes  the  results  of  four  rounds  of  directed  evolution  starting  from  lead 
variant  2D8,  and  screening  for  GB  and  GD  hydrolysis.  Screening  -2000  variants  in  round  1,  we 
identified  21  clones  that  were  improved  up  to  4-fold  with  GB,  and  up  to  5-fold  with  GD,  relative 
to  2D8,  in  crude  cell  lysates.  The  two  most  active  variants  with  GD,  5H8  and  PG11,  were 
purified  and  characterized.  We  found  improvements  of  21-  and  4-fold  in  the  catalytic  activity  of 
the  best  variant,  PG11,  with  the  two  toxic  isomers  of  GD  relative  to  the  starting  variant,  2D8. 
However,  its  activity  with  GF  was  reduced  2.3-fold  relative  to  2D8,  and  no  improvement  in  GB 
activity  occurred.  The  two  best  variants  shared  the  same  4  active-site  mutations:  Leu69Val, 
Hisl34Arg,  Phe222Ser  and  Thr332Ser,  but  differed  at  position  115:  5H8  acquired  a  Val  and 
PG11  an  Ala.  While  mutations  Phe222Ser  and  Thr332Ser  were  already  present  in  the  parent 
variant,  2D8,  the  mutations  in  positions  69,  115  and  134  were  selected  from  the  substitution 
library  of  round  1. 

For  the  2nd  round,  we  constructed  a  substitution  library  that  explored  various  substitutions  at 
position  222,  whilst  shuffling  the  10  most  active  clones  found  in  round  1.  The  resulting  library 
was  screened  on  GD  and  GB,  as  well  as  on  GF.  The  latter  was  included  to  eliminate  variants  that 
exhibit  reduced  activity  with  GF,  as  observed  with  variants  isolated  from  the  1st  round. 
Screening  of  -1000  clones  we  found  39  variants  that  were  improved  up  to  4-fold  towards  GB,  2- 
fold  towards  GD  and  4-fold  towards  GF  relative  to  the  best  variants  of  round  1.  Sequencing 
revealed  25  unique  variants  and  14  duplicated  variants.  Following  more  detailed  screens  (at 
different  agent  concentrations),  we  selected  the  6  most  active  variants,  purified  and  characterized 
them.  All  purified  variants  had  improved  catalytic  rates  up  to  7-  and  193-fold  towards  the  toxic 
isomers  of  GD,  and  up  to  5-fold  towards  GB,  relative  to  the  starting  variant,  2D8.  Their  activities 
towards  GF  were  improved  relative  to  variants  from  round  1,  and  they  became  similar  to  that  of 
2D8.  All  round  2  variants  contained  three  mutations:  Leu69Val,  Hisl34Arg  and  Thr332Ser,  and 
all  but  one  contained  also  the  Phe222Met  mutation.  The  greatest  variability  was  seen  at  position 
115,  where  either  Ala,  Leu  or  Val  were  observed.  While  no  attempt  to  introduce  random 
mutations  during  library  construction  was  made,  the  PCRs  applied  for  library  making  produced 
random  mutations  such  as  Phe64Leu  or  Asp309Gly,  that  were  retained  in  certain  improved 
variants.  By  round  2,  the  most  improved  variants  we  isolated  reached  the  targeted  catalytic 
efficiency  of  >10  [M’  min"  ]  for  GD  and  GF,  but  were  still  10-fold  lower  in  activity  towards 
GB.  The  best  round  2  variant,  VII-D11,  hydrolyzed  the  two  toxic  isomers  of  GD  with  equally 
high  rates  (2.9xl07  [M  1  min"1]. 

During  the  3  round,  two  libraries  were  constructed  and  screened  for  variants  with  increased 
activities  using  GD  and  GB.  The  libraries  were  made  by  shuffling  the  7  most  active  variants 
from  round  2,  together  with  oligonucleotides  encoding  for  specific  mutations  in  PONl’s  active 
site  and  flexible  active-site  loop.  The  construction  of  two  such  libraries  enabled  us  to  screen  a 
larger  diversity  of  mutations  in  separate  regions  of  the  enzyme.  Library  #1  introduced  mutations 
in  active-site  positions  that  had  not  yet  been  mutated  in  the  previous  rounds  and  in  second  shell 
positions  around  the  active  site,  as  well  as  ancestral  substitutions.  Library  #2  included  mutations 
in  the  active-site  loop.  Screening  400  variants  from  each  library  with  GD  and  GB  identified  a 
total  of  18  improved  clones.  All  selected  clones  were  improved  for  GB  but  only  67%  (12/18) 
were  also  improved  for  GD.  Most  of  the  improved  variants  originated  from  Library  #1  (73%, 
13/18),  and  the  improvements  for  GB  (1. 6-4-fold)  were  also  greater  than  for  GD  (1. 1-1.6  fold). 
Sequencing  indicated  that  55%  (10/18)  of  the  selected  variants  contained  at  least  one  ancestral 
mutation  (Leu55Ile,  or  Aspl36His)  and  28%  (5/18)  contained  two  ancestral  mutations.  In 
addition  to  these  ancestral  mutations,  the  active-site  mutation  Ile291Leu  became  abundant  in 
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selected  clones  (39%,  7/18).  Certain  mutations  introduced  in  earlier  rounds,  such  as  Hisl34Arg 
and  Phe222Met,  were  fixed  (/.<?.,  appeared  in  all  selected  clones)  and  others,  such  as  Leu69Val 
and  Hisll5Ala,  were  nearly  fixed  (90%,  16/18).  We  purified  and  characterized  the  four  most 
improved  variants  from  this  round.  The  catalytic  efficiencies  of  these  variants  were  improved  up 
to  11-  and  293-fold  for  the  two  toxic  isomers  of  GD,  and  up  to  17-fold  for  GB,  relative  to  2D8, 
our  starting  variant.  A  small  3-fold  improvement  was  found  for  GF  hydrolysis.  The  catalytic 
efficiency  of  our  most  improved  variant,  1-I-F11,  with  GB  (3.9xl06M  1  min"1)  was  improved  >3- 
fold  relative  to  the  best  round  2  variant,  but  was  still  2.5-fold  lower  than  our  goal. 

The  round  4  library  was  constructed  by  shuffling  the  1 8  most  improved  variants  from  round 
3  with  the  aim  of  recombining  beneficial  mutations  identified  in  separate  variants.  We  screened 
700  library  clones  using  GB  and  GD,  and  isolated  22  improved  variants  of  which  18  were 
unique.  Of  these,  66%  (12/18)  were  improved  mostly  for  GB,  and  33%  (6/18)  were  improved 
also  for  GD.  In  addition  to  the  mutations  fixed  in  earlier  rounds,  all  selected  Round  4  variants 
now  contained  the  mutation  Hisll5Ala,  and  most  variants  also  carried  Leu55Ile  (55%,  10/18) 
and/or  Ile291Leu  (83%,  15/18).  We  purified  9  variants,  and  identified  three  that  had  improved  up 
to  13-  and  340-fold  for  the  toxic  isomers  of  GD,  up  to  7-fold  for  GB,  and  up  to  9-fold  for  GF, 
with  respect  to  the  starting  point  variant  2D8.  Although  their  activities  with  GB  (1.7-3.2xl06  [M  1 
min"1])  were  3-6-fold  lower  than  our  desired  goal,  their  catalytic  efficiencies  with  GF  and  GD 
were  well  over  1x10  [M"  min"  ]. 

C.  Catalytic  activity  of  evolved  variants  with  GA 

The  in-vivo  toxicity  of  GA  is  >2-fold  lower  than  that  of  all  other  G-agents5,  thus  ranking  its 
toxicity  as  the  least  threatening  of  the  group.  However,  its  structure  (Figure  la)  is  significantly 
different  than  that  of  all  the  other  G-agents,  and  its  Rp  isomer  is  more  toxic  than  its  Sp  isomer. 
Therefore,  we  decided  to  examine  the  ability  of  wild-type  like  rePONl,  and  of  the  most 
improved  variants  from  each  round  of  evolution,  to  hydrolyze  GA.  Using  the  AChE  inhibition 
assay,  we  found  that  wt  rePONl  (G3C9)  is  3-30-fold  more  efficient  at  hydrolyzing  the  toxic 
isomer  of  GA  than  the  toxic  isomers  of  all  other  G-agents  (Table  7).  In  addition,  we  found  that 
our  most  improved  variants  from  rounds  3  and  4,  1-I-F11  and  IIG1,  were  4-5  fold  more  efficient 
than  rePONl  at  hydrolyzing  the  toxic  isomer  of  GA  (Table  7).  Thus,  in  spite  of  any  direct 
selection  pressure,  the  catalytic  activity  of  our  evolved  variants  with  GA  (1.7-2. 3  xlO6  M  1  min"1) 
had  increased  and  had  become  similar  to  their  activity  on  GB  (3. 2-3. 9  xlO6  M"1  min’1]). 

D.  Stereospecificity  of  the  evolved  variants 

Our  previous  efforts  to  evolve  PON1  for  hydrolysis  of  G-agents  made  use  of  their  coumarin 
analogues1.  In  these  compounds,  the  fluoride  leaving  group  is  replaced  by  a  coumarin  derivative 
(Figure  lb),  thus  permitting  fluorogenic  or  colorimetric  detection  of  enzymatic  activity  .  Due  to 
its  pKa  value  (3.1),  fluoride  is  a  better  leaving  group  than  3-cyano-7-hydroxy-4-methylcoumarin 
(pKa  =  6.51+0.04),  and  is  also  much  smaller.  Nevertheless,  the  activities  of  OP-hydrolyzing 
enzymes  on  the  coumarin  analogues  and  on  the  actual  G-type  nerve  agents  have  been  shown  to 
be  quite  similar  not  only  for  engineered  PON1  variants1,  but  also  for  other  OP-hydrolyzing 
enzymes,  such  as  PTE  .  In  this  project,  the  hydrolysis  of  G-agents  was  monitored  making  use  of 
the  AChE  assay  that  detects  only  hydrolysis  of  their  toxic  Sp  isomers.  To  examine  the  hydrolysis 
of  both  stereoisomers,  we  assayed  the  most  improved  variants  from  each  round  of  directed 
evolution  using  the  purified  Rp  and  Sp  isomers  of  the  coumarin  analogue  of  GF  (CMP-coumarin; 
Figure  lb).  The  results  indicate  that  the  increase  in  detoxification  rates  for  the  three  target  G- 
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agents  (GB,  GD  and  GF)  along  the  course  of  directed  evolution  was  accompanied  by  a  complete 
reversion  in  rePONl  stereoselectivity  in  favor  of  the  Sp  isomer  (Table  7).  The  catalytic  activity 
of  the  wild-type-like  rePONl  on  Rp-CMP-coumarin  is  >400-fold  higher  than  on  Sp-CMP- 
coumarin.  In  contrast,  the  catalytic  activity  of  the  best  variants  from  round  3  and  4  on  Sp-CMP- 
coumarin  was  found  to  be  >4000-fold  higher  than  their  activity  on  i?p-CMP-coumarin  (Table  8). 

We  were  able  to  attribute  the  increase  in  Sp/Rp  stereoselectivity  with  CMP-coumarin  to 
changes  both  in  substrate  binding  and  catalysis.  We  observed  a  gradual  decrease  in  Km  values  for 
Sp-CMP-coumarin  with  each  round  of  evolution  together  with  a  concomitant  increase  in  kcat 
values  (Table  8).  A  parallel,  gradual  increase  in  Km,  aaccompanied  by  a  decrease  in  kcat  values, 
was  observed  for  the  same  variants  with  Rp-CMP-coumarin.  The  change  in  the  enantiomeric 
ratio  (E)  differed  from  round  to  round,  with  the  most  dramatic  change  occurring  in  round  3 
variants.  A  similar  trend  was  observed  when  we  examined  the  hydrolysis  of  the  coumarin 
analogue  of  GD,  PMP-coumarin.  Unlike  other  G-agents,  GD  has  a  chiral  carbon  atom  in  addition 
to  its  chiral  phosphorus  center  (Figures  la  &  lb).  As  with  other  G-agents,  the  toxicity  of  GD  is 
determined  by  the  phosphorus  chirality,  and  GD's  two  Sp  isomers  (SpRc,  SpSc)  are  >  1000-fold 
stronger  inhibitors  of  AChE  than  its  two  Rp  isomers  ( RpSc,  RpRc)  ■  While  the  best  round  1 
variant,  PG11,  hydrolyzed  all  four  diasteroisomers  of  PMP-coumarin  at  measurable  rates, 
variants  from  rounds  2-4  hydrolyzed  only  50%  of  the  racemic  mixture.  To  confirm  that  the 
hydrolyzed  fraction  corresponds  to  the  Sp  isomer  pair,  we  made  use  of  the  evolved  rePONl 
variant  3B3,  that  had  previously  been  shown  to  hydrolyze  almost  exclusively  the  Rp  isomers  of 
various  coumarin  OPs,  including  (9-ethyl  methylphosphonyl-coumarin,  (9-n-  or  (9-i-propyl 
methylphosphonyl-coumarin,  and  (9-cyclohexyl  methylphosphonyl-coumarin  (CMP- 
coumarin)  .  Here,  we  observed  that  3B3  was  able  to  hydrolyze  only  50%  of  a  racemic  mixture 
of  PMP-coumarin,  and  that  the  remaining  50%  was  readily  hydrolyzed  by  evolved  variants  from 
rounds  2-4  (not  shown).  We  thus  conclude  that  in  response  to  the  continuous  application  of  a 
screen  for  protection  of  AChE,  which  identifies  enzymes  that  target  the  Sp  isomers  of  the  applied 
G-agents,  our  evolved  variants  became  highly  specialized  for  Sp  isomer  hydrolysis 

E.  Generation  of  PON1  variants  specific  for  V-type  nerve  agents 

The  evolution  of  V-type  hydrolyzing  variants  has  lagged  behind  the  development  of  G-type 
hydrolyzing  variants  primarily  due  to  the  lack  of  suitable  V-agent  surrogates.  Variants  that  failed 
to  hydrolyze  the  O,  (9-diethylphosphoryl  analog  of  V-type  nerve  agents,  amiton 
(C2H50)2P(0)SCH2CH2N(Et)2  (Fig.  4),  were  subsequently  found  to  exhibit  low,  yet  detectable, 
activity  on  the  (9-alkyl  methylphosphonothiolates  (e.g..  VM,  VX,  RVX).  Variants  belonging  to 
earlier  generations  (e.g.,  wt-G3C9,  3B3,  1A4,  0C9,  2D8),  hydrolyzed  VM  at  2-5  mU/mg 
compared  to  40  mU/mg  displayed  by  bacterial  wt  PTE.  Amiton  is  hydrolyzed  by  PTE  (50 
mU/mg)  but  not  by  the  rePONl  variants  mentioned  above.  Thus,  VM,  rather  than  amiton,  can 
serve  as  a  suitable  VX  analog  for  identification  of  rePONl  variants  capable  of  hydrolyzing  P-S- 
alkyl  containing  nerve  agents.  In  order  to  develop  an  efficient  screening  procedure  for 
identification  of  variants  capable  of  efficiently  degrading  V  agents,  we  again  utilized  the 
enzymic  approach  with  hAChE  added  to  the  lysates.  For  this  purpose,  we  developed  an  in  situ 
non-hazardous  synthetic  procedure  for  the  generation  of  two  pairs  of  methylphosphonothiolates 
that  differ  in  the  size  of  the  (9-alkyl  moiety,  namely  VX  and  its  homolog  VM,  and  RVX  and  its 
homolog  RVX-isoPr  (Fig.  4).  Our  first  (Round  1)  VX  library  was  made  by  screening  lysates  of 
the  R4  G  agent  library  (Table  6b)  with  VX,  utilizing  the  ISA  protocol  based  on  inhibition  of 
hAChE.  We  isolated  and  characterized  3  improved  variants  from  that  round,  which  displayed 
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>200-fold  enhanced  hydrolysis  of  VX  relative  to  wt-rePONl  (Table  9).  We  then  used  VX  to 
assay  the  catalytic  efficiency  of  variants  isolated  and  purified  from  all  previous  rounds  of  G- 
agent  evolution.  Although  variants  from  earlier  rounds  of  G-type  evolution  were  less  active  on 
VX,  they  had  different  mutational  compositions.  We  therefore  combined  the  4  variants  from 
earlier  rounds  with  the  highest  activity  towards  VX  with  the  3  most  active  R4  variants  referred  to 
above,  which  were  denoted  VX-R1,  and  shuffled  them  together  to  generate  a  Round  2  library  of 
VX  mutants  (VX-R2  library). 

The  round  2  library  was  screened  again  against  VX  and  RVX,  and  it  was  possible  to  isolate 
and  characterize  4  variants  that  showed  significantly  improved  VX  hydrolysis.  The  best  VXR2 
variant  was  S-2-C8  (Table  9).  In  conclusion,  efforts  towards  achieving  this  milestone  have 
resulted,  so  far,  in  at  least  one  variant  capable  of  hydrolyzing  the  toxic  isomer  of  VX  >5000-, 
>130-  and  >25-fold  better  than  the  wt-G3C9,  4E9  and  the  VXG1  variants,  respectively  (Table 
9).  Thus,  S-2-C8  hydrolyzed  the  toxic  isomer  of  VX  at  l.OxlO4  M  1  min"1,  a  value  similar  to  that 
for  a  bacterial  PTE  that  we  had  used  to  verify  the  interception  protocol  with  V  agents  (Table  9). 
To  gain  further  insight  into  the  stereo-preference  of  S-2-C8  when  reacting  with  the  racemic  VX, 
release  of  the  leaving  group,  HSCH2CH2N(isoPr)2,  was  monitored  by  use  of  Ellman’s  reagent, 
DTNB.  Fig.  5  clearly  shows  that  the  two  optical  isomers  of  VX  are  hydrolyzed  equally  well  by 
S-2-C8.  A  similar  kinetic  time  course,  albeit  slower,  was  observed  with  VM  (not  shown).  RVX 
and  RVX-isPr  were  found,  however,  to  be  poor  substrates,  suggesting  that  although  the  binding 
pocket  in  this  variant  accommodates  the  (9-ethyl  moiety  well,  it  needs  to  be  enlarged  to  permit 
enhanced  hydrolysis  of  the  (9-isobutyl-containing  V  agents. 

To  further  improve  the  VX  hydrolysis  capacity  of  PON1  we  screened  a  round  3  VX  library 
VX  (VXR3),  which  was  made  by  shuffling  the  10  best  clones  of  round  2  with  oligos  introducing 
34  new  site-specific  mutations.  We  identified  19  improved  variants  from  -1000  random  clones. 
Then  we  purified  and  characterized  4  of  them.  The  best  variant,  1-1-H3,  was  improved  by  3.5- 
fold  over  the  best  variant  of  the  previous  round,  S-2-C8,  and  had  kcat/Km  -35,000  [M  'min'1] 
towards  the  toxic  isomer  of  VX.  Notably,  when  the  rate  of  release  of  the  leaving  group, 
HSCH2CH2N(iPr)2,  was  determined  using  DTNB,  kcat/km  for  1-1-H3  was  slightly  lower, 
~2.4xl04  IVr'min"1,  suggesting  that  the  two  isomers  of  VX  are  hydrolyzed  at  different,  albeit  low, 
rates.  When  detoxification  of  RVX  was  followed  by  the  enzymic  protocol,  1-1-H3  displayed 
poor  activity,  and  could  not  be  determined;  however,  hydrolysis  of  the  second  half,  i.e.,  of  the 
less  toxic  isomer,  as  evaluated  using  DTNB,  was  -700  M  1  min"1.  VM,  the  N,N, -diethyl  analog  of 
VX  was  hydrolyzed  at  approximately  lxlO4  M-lmin-1,  thus  indicating  that  the  steric  constraint 
presented  by  the  isobutyl  moiety  of  RVX  does  not  permit  proper  accommodation  in  the  catalytic 
site.  Results  from  screening  the  VXR3  library  were  taken  into  consideration  in  building  the 
VXR4  library. 

We  made  three  libraries  of  variants  by  shuffling  the  10  best  clones  from  round  VXR3  and 
introducing  mutations  at  4  new  positions.  The  three  libraries  differed  in  the  number  of  new 
mutations  introduced  per  variant:  1,  1.5  and  2.5  on  average.  Novel  mutations  were  observed  in  2 
positions.  The  preliminary  results  from  the  enzymic  protocol  processed  in  well  lysates  show  a 
-10-fold  increase  in  proficiency  towards  VX  when  compared  to  VXR3  library.  The  most  active 
VXR4  variant  displayed  kcat/Km=1.5xl05  M  'min"1,  >75,000-fold  greater  than  the  wt  G3C9 
enzyme. 
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F.  Ex-vivo  experiments  in  human  whole  blood2 

Ex  vivo  protocol:  In  order  to  evaluate  the  stability  and  proficiency  of  rePONl  variants  in 
human  blood  we  developed  an  ex-vivo  protocol  to  screen  the  potential  antidotal  capacity  of 
rePONl  variants  via  their  protection  of  endogenous  acetylcholinesterase  (AChE)  and 
butyrylcholinesterase  (BChE)  in  human  whole  blood  obtained  from  a  local  blood  bank.  This 
protocol  permitted  us  to  address  the  relationship  between  blood  rePONl  concentrations,  together 
with  their  kinetic  parameters,  and  the  level  of  protection  conferred  on  the  ChEs  in  human  blood 
following  a  challenge  with  GF.  The  GF  was  generated  in  situ  at  a  non-hazardous  concentration. 
In  addition,  long-term  incubation  of  the  tested  rePONl  variants  with  whole  blood,  prior  to 
spiking  with  GF,  provided  information  concerning  the  effect  of  their  association  with  blood 
components  on  their  stability.  For  example,  ex-vivo  treatment  of  human  whole  blood  with  0.45 
pM  VIID2  (an  R4  variant)  resulted  in  retention  of  42%  of  the  total  ChE  activity  following 
spiking  with  0.08  pM  GF,  compared  to  only  5%  in  unprotected  blood.  These  experimental  data 
are  in  good  agreement  with  the  %  residual  activities  calculated  on  the  basis  of  the  rate  constants 
for  inhibition  of  human  AChE  and  BChE  by  GF,  the  concentration  of  the  rePONl  variant,  and  its 
kcat/Kra  value.  The  validation  protocol  developed  provides  a  rapid  and  reliable  ex-vivo  screening 
tool  for  selection  of  rePON  1  bioscavenger  candidates  suitable  for  use  in  the  protection  of  humans 
against  OP  intoxication.  The  results  also  provide  input  for  a  revised  mathematical  model  for 
estimating  the  efficacious  dose  required  in  humans  of  a  given  rePONl  variant. 

The  use  of  evolved  PON1  variants  as  prophylactic  drugs  requires  that  they  maintain  their 
catalytic  activity  in  the  complex  chemical  and  biological  environment  of  human  blood  over  long 
periods  of  time.  Previously,  we  showed  that  evolved  variants  4E9  and  VIID2  were  able  to  protect 
human  blood  cholinesterases  (ChE)  from  inhibition  by  GF  using  an  ex-vivo  assay34.  Here  we 
examined  the  duration  of  blood  ChE  protection  from  GF  using  representative  variants  from  each 
round  of  selection.  In  order  to  do  so,  we  added  purified  variants  PG11,  VIID11,  1-I-F11  and 
VIID2  (from  rounds  1-4  respectively),  to  human  blood  samples,  followed  by  incubation  at  37°C 
for  24  h.  Samples  of  blood-PONl  variant  mixtures  were  then  spiked  with  in  situ  generated  GF  at 
different  time  intervals,  and  their  residual  ChE  activity  was  then  measured.  The  ability  of  each 
variant  to  protect  blood  ChE  from  inhibition  by  GF  at  each  time  point  was  then  calculated 
relative  to  its  initial  protection  activity  obtained  after  only  5  min  of  incubation  in  blood  (Figure 
6,  upper  panel).  While  round  1  variant  PG11  could  not  provide  measurable  ChE  protection  at 
the  PON  concentration  range  used  (0.5-2. 1  pM)  due  to  its  low  catalytic  activity,  variants  VIID11, 
1-I-F11  and  VIID2  were  able  to  maintain  80±2  %,  95±12  %  and  95±2  %  of  their  initial  ChE 
protection  activity  in  human  blood  following  24  h  incubation  at  37  °C.  The  ChE  activity  of  blood 
samples  not  containing  PON  variants  exhibited  no  detectable  changes  following  24  h  of 
incubation  at  37 °C  (data  not  shown). 

We  then  examined  the  stability  of  our  evolved  variants  under  the  same  conditions  by 
assaying  their  ability  to  hydrolyze  a  chromogenic  analogue  of  GF,  CMP-coumarin,  at  various 
time  intervals  (Figure  6.  Lower  panel).  This  assay  provides  an  assessment  of  the  variants 
stability  by  directly  reporting  on  their  catalytic  activity.  Samples  of  blood-PONl  mixtures 
incubated  at  37  °C  were  diluted  1000-fold  into  buffered  solutions  containing  CMP-coumarin,  and 
the  rate  of  substrate  hydrolysis  (i.e.,  release  of  the  coumarin  leaving  group)  was  monitored  for  5 
min  at  400  nM.  The  results  were  plotted  as  percent  of  the  initial  activity  at  t  =  0  (Figure  6,  lower 
panel).  As  with  the  previous  assay,  variants  VIID11,  1-I-F11  and  VIID2  maintained:  77±12%, 
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"  For  experimental  details  see  Ashani  et  al  (201 1)  ref  #34. 
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95±6  %  and  92±5%  of  their  initial  activities  following  a  24  h  incubation  period  in  human  blood 
samples  ( ex-vivo )  at  37°C.  The  stability  of  variant  PG11  (76±5  %)  was  similar  to  that  of  1-I-F11. 
Surprisingly,  the  activity  of  all  variants  after  1  h  incubation  increased  4-12%  relative  to  their 
activity  after  5  min.  The  spontaneous  hydrolysis  of  CMP-coumarin  by  blood  samples  devoid  of 
any  PON1  variant  was  found  to  be  negligible  (data  not  shown).  The  results  from  the  ex-vivo 
experiments  were  similar  to  the  stability  of  the  4  best  variants  observed  in  Tris  buffer,  pH  7.4 
(37°C)  (not  shown).  Thus,  taken  together,  these  results  indicate  that  the  ex  vivo  interaction  of  the 
evolved  PON  1  variants  with  human  blood  constituents  did  not  affect  their  proficiency  at  least  for 
24  h  at  37°C. 

Key  Research  Accomplishments 

•  Development  of  a  screen  for  PON1  variants  with  increased  OP  activity,  based  on  sorting 
emulsion  droplet  by  FACS. 

•  Development  of  a  novel  neutral  drift  PON  1  library 

•  Synthesis  of  amiton,  a  P-S  model  compound  for  VX,  and  safe  and  non-hazardous  procedures 
for  in  situ  generation  of  G-  and  V-type  nerve  agents 

•  Protocol  for  stereo-specific  synthesis  of  a  chiral  (9-alkyl  methylphosphonate  with  a  coumarin 
fluorescent  leaving  group 

•  Interception  protocol  (ISA)  for  direct  screening  of  crude  lysates  with  G-type  nerve  agents 

•  Isolation,  overexpression  and  purification  of  several  variants  with  greatly  enhanced  hydrolysis 
towards  GA,  GB,  GD.  GF,  (kcat/Km  of  2-5xl07  M  'min"1)  and  towards  VX  (kcat/Km=1.5xl05 
Mimin'1 ). 

•  Libraries  of  variants  capable  of  hydrolyzing  G-  and  V-type  nerve  agents 

•  3D  structures  of  wt  rePONl  and  of  5  additional  catalytically  active  PON1  variants. 

•  A  protocol  for  ex-vivo  estimation  of  kcat/Km  values  of  rePON  1  variants  in  human  blood 

•  Several  variants  were  over-expressed,  purified  and  transferred  along  with  their  plasmid  to 
USAMRICD  (Aberdeen  Proving  Ground,  MD)  and  to  the  Department  of  Biochemistry,  Ohio 
State  University,  for  further  in  vivo  evaluation  of  protection  against  G-type  nerve  agent 
intoxication  (including  variants:  2C3,  3D8,  8C8,  1G3,  3B3,  2H4,  0C9,  2D8  and  1A4,  4E9, 
MG2-1A4,  VIID2). 
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Conclusions 

Directed  evolution,  combining  random  and  designed  mutations  based  on  3D  structures, 
resulted  in  the  identification,  expression  and  purification  of  mutants  of  a  recombinant 
mammalian  PON1  (rePONl)  with  catalytic  activity  that  approaches  the  theoretical  estimated 
minimal  requirement  for  efficient  prophylactic  protection  at  reasonable  protein  doses.  The 
kcat/Km  values  for  these  variants,  when  reacting  with  G-type  CWNAs,  are  as  high  as  2-5xl07  M 
'min'1,  and  it  is  envisaged  that  this  proficiency  will  be  further  increased  by  additional  rounds  of 
directed  evolution,  and  thus  further  decrease  the  loading  dose  of  the  protecting  enzyme.  The 
most  advanced  mutants  were  evolved  to  completely  reverse  the  stereo-preference  of  PON1  from 
the  non-toxic  to  the  toxic  isomers  of  racemic  nerve  agents.  PON1  variants  evolved  by  selection 
with  VX  are  capable  of  hydrolyzing  the  toxic  isomer  of  VX  >75, 000-fold  faster  than  the  wt 
enzyme,  with  kcat/Kra=1.5xl05  M  'min  '.  Additional  rounds  will  be  required  to  obtain  the 
efficiency  required  to  qualify  PON1  as  antidotes  against  VX  (kcat/Km  >lxl07  M  'min  ').  These 
results,  taken  together  with  preliminary  in  vivo  evaluation  (performed  at  ICD,  Aberdeen  Proving 
Ground  MD),  suggest  that  one  or  two  PON1  variants  might  serve  as  broad- spectrum  catalytic 
bioscavengers  to  protect  against  multiple  LD50  doses  of  all  realistic  CWNA  threat  agents  at  low 
protein  doses.  This  technology  is  further  envisaged  to  provide  rapid  discovery  of  pretreatment 
and  post-challenge  therapeutic  drugs  not  only  against  existing  threats  but  also  to  counteract 
emerging  OP-based  CWNA  and  commercial  pesticides,  and  thus  shorten  the  time  from 
emergence  of  a  threat  to  identification  of  potential  antidotes  measures  to  a  few  days  or  weeks. 
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Tables 


Table  1.  Activity  of  P0N1  mutants  with  both  isomers  of  CMP-coumarin  and  with 
Sp  -IMP-coumarin 


Variants 

Sp-CMP 

coumarin  a 

(kCat/Km) 

Mimin'1 

Sp-IMP 

coumarin a 

(kcat  /Km) 

Mimin'1 

Rp-CMP 

coumarin  a 
*106  Apparent 

(kcat/Kn) 

Mimin'1 

Non- 

synonymous 

mutations0 

Wild-type-like 

<200 

n.d. 

0.08±0.0034 

- 

rePONI  -G3C9 

(1)b 

(1) 

H115W 

331 ±39 

1 983±30 

0.45 

H115W 

(>1.6) 

(1) 

(6)t 

V346A 

885±76 

3633±1 26 

0.2 

V346A 

(>4.4)  t 

(1  -8)t 

(2.5)  t 

HI  15W+V346A 

81 3±79 

1 01 40±7 

0.4 

H115W, 

(>4-1)T 

(5)t 

(5)  T 

V346A 

a.  For  each  variant,  enzymatic  activities  (kcat  /Km)  were  measured  with  purified  proteins  and 
denoted  are  the  average  ±  standard  deviation  values  obtained  from  3  independent  repeats.  The 
values  without  standard  deviations  had  s.d  <20%  of  their  values.  Values  in  parentheses 
denoted  the  fold  increase  and  decrease  as  compared  to  wt-like  rePONI  for  either  isomer  of 
CMP  and  as  compared  to  HI  15W  for  SP-IMP.  n.d.  denotes  non-detectable  activity. 

b.  The  catalytic  efficiency  was  estimated  as  described  in  Ref#1 

c.  Denoted  in  bold  are  mutations  in  active-site  residues. 
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Table  2:  Representative  variants  along  the  directed  evolution  process 


Variant  3 

Mutations  b 

Sp-CMP-Coumarin  c 

A?P-CMP-Coumarin 
Apparent  kcat/  KM  c 
(pM^min'1) 

Sp-CMP-F 

Apparent 

kCat/KMc  (pM^min'1) 

kCat 

(min’1) 

Km 

(MM) 

kcat/  Km 

(|jM'1min‘1) 

rePON1G3C9 

(Wild-type-like) 

nd 

nd 

<0.0002 

(1)d 

0.08±0.0034  (1) 

0.1 3±0.03 

(D 

3B3 

N41D,  SHOP, 

L240S,  H243R, 

F264L,  N324D, 

T332A 

nd 

nd 

<0.0002 

(1)d 

20±1 .7 

(250)  t 

-0.0001 

H115W- 

V346A 

H115W,  V346A 

nd 

nd 

0.0008 

(>4)t 

0.4 

(5)  T 

0.02 

±0.003 

3A7 

V97A,  H115W, 

P135A,  F222S, 

M289I 

nd 

nd 

0.0027 

(>13.5)  t 

0.16 

(2)1 

0.008 

±0.0035 

8C8 

L69S,  V97A, 

H115W,  P135A, 

F222S 

11.6 

±0.18 

124.5 

±5.8 

0.093  ±0.003 

(>465)  t 

0.0035 

(23)  i 

0.2 

(1.5) 

2D8 

L69G,  H115W, 

H134R,  F222S, 

T332S 

268 

±1.6 

76.3 

±3.2 

3.52 

±0.13 

(>17600)  t 

0.465 

(5.8)  1 

14.3 

(110) 

3D8 

L69G,  H115W, 

H134R,  Ml  96V 

F222S,  T332S 

295 

±1.62 

25.4 

±0.5 

11.6 

±0.23 

(>58000)  t 

nde 

3.3 

(25) 

4E9 

L69G,  S111T, 

H115W,  H134R, 

F222S,  T332S 

513 

23 

22.3 

(>111500) 

T 

nde 

16.8 

(129) 

a.  Annotation  of  variants:  The  first  digit  relates  to  the  plate  number,  and  the  following  letter-digit  to  its 


location  within  the  plate.  For  example,  variant  3A7  =  plate  #3,  well  A7;  n.d.,  not  detectable. 

b.  Denoted  in  bold  are  mutations  in  active-site  residues. 

c.  Enzymatic  parameters  were  measured  with  purified  proteins  and  comprise  the  average  obtained 
from  the  3  independent  repeats.  Error  ranges  represent  the  standard  deviations  observed  between 
measurements.  Values  in  parentheses  describe  the  fold-change  compare  to  the  starting  point, 
rePON-G3C9.  The  kinetic  parameters  for  Sp-CMP-coumarin  were  spectrophotometrically  measured 
with  pure  substrate  samples.  Parameters  for  fiP-CMP-coumarin  were  determined  with  the  racemate, 
and  for  CMP-F  with  the  in  situ  prepared  substrate  and  an  AChE  inhibition  assay. 

d.  The  catalytic  efficiency  was  estimated  as  described  in  Ref.1. 
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Table  3.  Improved  3rd  round  variants  from  libraries  derived  from  rePON1-H115W- 
V346A. 


Variants  a 

Fold 

improvement 
with  Sp-IMP- 
coumarin b 

Fold 

improvement 
with  Sp-CMP- 
coumarin b 

Non-Synonymous 
mutations  c 

4D2 

0.3x 

lOx 

L69S,  H115W,  P135A,  F222S 

8C8 

0.4x 

1 3x 

L69S,  V97A,  H115W,  P135A, 
F222S 

6C5 

2x 

3x 

V97A,  H115W,  P135A, 

F222S,  Ml 96V,  M289I 

1A8 

1.3x 

1.5x 

L4P,  V97A,  H115W,  P135A, 
F222S,  Ml 96V 

1E3 

0.2x 

1.4x 

A6T,  V97A,  H115W,  P135A, 
F222S,  D212N,  M289I 

7G10 

0.5x 

1.3x 

LI  OS,  H115W,  F222S,  M289I, 
V346A 

8H4 

0.7x 

0.8x 

V97A,  H115W,  P135A, 

F222S 

8H3 

0.8x 

0.8x 

V97A,  H115W,  P135A, 

F222S,  1237V,  L262F 

1G1 

1.2x 

1.2x 

V97A,  H115W,  F222S,  M289I 

a.  The  annotation  of  the  variants:  The  first  letter  relates  to  the  plate  number,  and  the  letter-digit  to 
the  location  of  the  clone  within  this  plate.  For  example,  variant  4D2=  plate  #  4,  well  D2. 

b.  Shown  are  all  variants  that  exhibited  higher  Sp  -IMP-coumarin  and  Sp-CMP-coumarin  activities  in 
crude  lysates  relative  to  the  H115W  +  V346A  PON1  mutant.  For  each  variant,  enzymatic 
activities  were  measured  in  crude  lysates,  and  denoted  are  the  average  values  of  fold 
improvement  obtained  from  3  independent  repeats.  The  values  had  s.d.<20%  of  their  value. 

c.  Non-synonymous  mutations  observed  in  each  variant.  Mutations  in  active-site  residues  are  noted 
in  red. 
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Table  4.  Improved  variants  from  shuffling  of  the  targeted  substitutions  library 
(6th  round  variants)). 


CMP(Sp) 

75 

CMPD 

(Ftp) 

CMPD 

(racemic) 

IMPD 

(Sp) 

EMPD 

(Sp) 

EMPD 

(Sp) 

Variant 

kcat 

Km 

kcat  /Km 

kcat  /Km 

Non-synonymous 

a 

min'1 

|JM 

IjM^min'1 

|jM' 

^min'1 

mutations  c 

8C8 

11.6 

124.5 

0.093 

0.0035 

0.03(1) 

0.022 

0.08 

0.2 

L69S,  V97A,  H115W, 

±0.18 

±5.8 

±0.003 

(1) 

(1) 

(1) 

(1) 

P135A,  F222S 

(1) 

2C3 

149.5 

212.7 

0.7 

0.088 

0.62 

0.079 

0.88 

0.16 

L69G,  H115W,  H134R, 

±0.87 

±5.5 

±0.01 

(25)  T 

(21)  T 

(3.6) 

(11)  T 

(1) 

F222S,  K233E 

(7.5)  T 

T 

5H5 

126 

102 

1.25 

0.2088( 

2.43 

0.76 

3.8 

0.8 

LI  OS,  F28Y,  L69G, 

±2 

±3.8 

±0.05 

60)  T 

(81)  T 

(35)  T 

(48)  T 

(4)  T 

H115W,  H134R,  F222S, 

(13.4)  T 

T332S 

0C9 

185 

65 

2.85 

0.296 

3.56 

1.04 

6.67 

2.6 

L14M,  L69G,  SUIT, 

±2.5 

±3.6 

±0.1 

(85)  T 

(119)  T 

(47)  T 

(83)  T 

(13)  T 

H115W,  H134R,  F222S, 

(31)  T 

T332S 

2D8 

268 

76.3 

3.52 

0.465 

3.04 

0.85 

6.8 

0.98 

L69G,  H115W,  H134R, 

±1.6 

±3.2 

±0.13 

(133)  T 

(101)  T 

(39)  t 

(85)  T 

(5)  T 

F222S,  T332S 

(38)  T 

1A4 

185 

51 

3.63 

0.39 

3.11 

0.85 

7.15 

2.2 

A6E,  L69G,  H115W, 

±1.5 

±1.7 

±0.1 

(1 1 1)T 

(104)  t 

(39)  t 

(89)  T 

(11)  T 

H134R,  F222S,  K233E, 

(39)  T 

T332S,  T326S 

3D8 

295 

25.4 

11.6 

nda 

4.73 

4.6 

12 

6.2 

L69G,  H115W,  H134R, 

±1.62 

±0.5 

±0.23 

(158)t 

(209) 

(150) 

(31)  T 

Ml 96V,  F222S,  T332S 

(125)  t 

T 

T 

a.  The  annotation  of  variants:  The  first  letter  relates  to  the  plate  number,  and  the  letter-digit  to  the 
location  of  the  clone  within  this  plate.  For  example,  variant  3B3  =  plate  #  3,  well  B3. 


b.  For  each  variant,  enzymatic  activities  (kcat/Km)  were  measured  with  purified  proteins,  and  denoted 
are  the  average  ±  standard  deviation  values  obtained  from  the  3  independent  repeats.  The 
individual  values  exhibited  s.d  <20%.  Values  in  parentheses  denoted  the  fold  increase  and 
decrease  as  compare  to  8C8,  the  best  variant  of  the  previous  round. 

c.  Denoted  in  bold  are  mutations  in  active-site  residues. 

d.  Variant  exhibited  single-phase  kinetics  of  product  release  when  reacted  with  racemic  CMP- 
coumarin,  suggesting  that  the  rates  of  hydrolysis  for  RP-  and  Sp-coumarin  are  similar. 


Table  5:  Activities  of  selected  rePONI  variants  on  Sp-CMP-coumarin  and  its 
fluoridated  product  CMP-F  (cyclosarin).a,b,c’d 


mutant 

(Sp)CMP- 

coumarin 

CMP- 

fluoridate 

fluoridate/ 

coumarin 

8C8 

0.09 

0.2 

2.2 

3D8 

11.6 

3.3 

0.3 

0C9 

2.8 

11.1 

3.9 

2D8 

3.5 

14.3 

4.1 

1A4 

3.6 

11.3 

3.1 

2C3 

0.7 

0.47 

0.7 

a.  The  figures  shown  are  values  of  kcat/KMx1 0b  (Mimin'1) 

b.  Data  for  OP-coumarin  are  based  on  release  of  the  chromophore  monitored  at  400  nm. 

c.  The  kcat/Km  values  for  the  fluoridates  were  determined  by  monitoring  the  rate  of  loss  of  anti-AChE 
potency  of  the  in  s/'fu-generated  compound,  assuming  Km»[P-F],  Calculations  are  based  on  a 
single  enzyme  concentration  selected  to  fit  the  dynamic  range  for  determination  of  the  apparent 
kobs  of  loss  of  anti-AChE  potency. 

d.  The  coumarin  leaving  group  was  replaced  by  fluoride  in  racemic  CMP-coumarin  to  yield  the 
racemic  fluoridate  of  CMP  (CMP-F).  It  should  be  noted  that  the  data  for  the  hydrolysis  of  CMP-F 
can  be  attributed  mostly  to  its  toxic  (SP)  isomer. 
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Table  6a.  Sequences  and  fold  improvement  of  selected  variants  from  round  1 


A.  Sequences  and  fold  improvement  with  sarin. 


3.9  3.5  3.5  2  1.7  1.6  1.4  1.4  1.3  1.2  1.0  1.0  0.9  0.8  0 


rePONI 

Amino 

Acid 

2D8 

9C3 

4F3 

5B2 

PIH12 

PD4 

PG3 

PA2 

IG6 

3D6 

IB3 

5H8 

9H3 

JD4 

PG11 

E2B 

69 

L 

G 

V 

V 

V 

L 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V  i 

70 

K 

K 

K 

K 

K 

N 

K 

K 

Q 

K 

K 

Q 

K 

K 

K 

K 

K 

71 

Y 

Y 

Y 

Y 

i 

i 

Y 

i 

Y 

Y 

Y 

i 

Y 

Y 

v  1 

82 

S 

s 

S 

s 

s 

s 

s 

s 

S 

s 

S 

S 

S 

s 

S 

S 

p 

115 

H 

w 

W 

L 

L 

L 

V 

V 

L 

V 

L 

L 

V 

w 

L 

A 

C 

134 

H 

R 

R 

R 

R 

R 

R 

R 

N 

R 

R 

N 

R 

R 

N 

R 

R 

196 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

L 

M 

M 

222 

F 

s 

s 

S 

s 

5 

5 

s 

S 

s 

s 

S 

S 

S 

5 

S  1 

240 

L 

L 

L 

L 

L 

L 

I 

I 

L 

I 

L 

L 

L 

L 

L 

L 

V 

292 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

L 

F 

F 

F 

F 

F 

F 

332 

T 

s 

5 

s 

s 

s 

s 

s 

S 

5 

5 

s 

s 

5 

*  First  row  -  Fold  Improvement  in  activity  of  each  variant  relative  to  control  variant  2D8  as  measured  in  cleared  cell 
lysates.  Numbers  represent  an  average  of  2  measurements  with  S.D.  <10%  of  value. 

**  Second  row  -  variant  names. 


B.  Sequences  and  fold  improvement  with  soman. 


5.2  4.9  3.7  3.3  3.0  2.5  2.5  2.5  2.5  2.4  2.4  2.3  2.3  2.2  2.2  1.7  1.4  1.2  1.0 


Amino 

rePONI  Acid 

2D8 

PIH12 

PG11 

5H8 

PG3 

1G6 

PD4 

E2B 

3D  6 

4G8 

4F3 

PA2 

5B2 

6C8 

2D1 

PF4 

IB3 

5A9 

JD4 

7A6 

69 

L 

G 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

V 

L 

70 

K 

K 

K> 

K 

K 

K 

K 

K 

K 

K 

K 

K 

9 

K 

5 

K 

K 

9 

K 

K 

* 

71 

Y 

Y 

Y 

Y 

Y 

i 

i 

i 

Y 

Y 

Y 

(BUI 

Y 

Y 

Y 

Y 

Y 

D 

Y 

Y 

82 

s 

s 

s 

s 

S 

s 

s 

s 

P 

s 

s 

s 

s 

s 

S 

P 

s 

S 

S 

s 

115 

H 

A 

V 

V 

V 

V 

c 

L 

V 

L 

L 

L 

V 

A 

c 

1. 

V 

L 

■■ 

134 

H 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

|  R 

N 

R 

R 

R 

R 

N 

R 

N 

R  zi 

196 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

L 

[  F  | 

M 

M 

M 

M 

M 

M 

M 

l 

M 

222 

F 

■■ 

s 
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s 

s 

s 

s 

s 

S 

S 

s 

s 

s 

s 

s 

s 

s 

S 

WM 

227 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

s 

L 

L 

240 

L 

L 

L 

L 

L 

I 

i 

i 

V 

L 

L 

L 

L 

L 

L 

V 

V 

L 

L 

L 

i 

292 

F 

F 

F 

F 

F 

F 

F 

F 

F 

L 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

309 

D 

0 

D 

D 

D 

D 

D 

D 

D 

D 

G 

D 

D 

D 

D 

D 

D 

D 

D 

D 

D 

332 

T 

5 

S 

s 

s 

5 

S 

5 

S 

S 

: 

S 

5 

S 

S 

S 

s 

DL2 

*  First  row  -  Fold  Improvement  in  activity  of  each  variant  relative  to  control  variant  2D8  as  measured  in  cleared  cell 
lysates.  Numbers  represent  an  average  of  2  measurements  with  S.D.  <10%  of  value. 

**  Second  row  -  variant  names. 


-25- 


Table  6b:  Catalytic  activities  of  improved  variants  from  each  round  of  evolution 


Variant 

Mutations  relative  to 
rePONI  a 

Round 

Catalytic  efficiency  b,d  (kCat/KM)  xIO7  Mimin'1 

GDC 

Fast  Slow 

GF 

GB 

re  PON  1 

- 

- 

0.0055 

±0.0017 

0.0015 

±0.0006 

0.01 3e 
±0.003 

0.008 

±0.001 

8C8 

Leu69Ser,  Val97Ala, 
Hisl  1 5T rp,  Pro135Ala, 
Phe222Ser 

0 

0.0028 

±0.0006f 

(0.5) 

0.0014 

±0.0004f 

(0.9) 

0.01 5e 
±0.005 
(1.2) 

0.0034 

±0.001 

(0.4) 

0C9 

Leu14Met,  Leu69Gly, 
Seri  1 1  Thr,  Hisl  1 5T rp, 
His134Arg, 
Phe222Ser, 
Thr332Ser 

0 

0.34 

±0.03f 

(62) 

0.026 

±0.0006f 

(17) 

1.11° 

±0.3 

(85) 

0.028 

±0.009 

(4) 

1A4 

Ala6Glu,  Leu69Gly, 
Hisl  1 5T rp,  His134Arg, 
Phe222Ser, 
Lys233Glu, 

Thr326Ser,  Thr332Ser 

0 

0.41 

±0.07f 

(75) 

0.033 

±0.004f 

(22) 

1.13e 

±0.3 

(87) 

0.023 

±0.003 

(3) 

4E9 

Leu69Gly,  Seri  1 1  Thr, 
Hisl  1 5T rp,  His134Arg, 
Phe222Ser, 
Thr332Ser 

0 

0.74 

±0.3 

(135) 

0.056 

±0.01 

(37) 

1.75e 

±0.3 

(135) 

0.033 

±0.004 

(4) 

5H8 

Leu69Val,  His115Val, 

His134Arg, 

Phe222Ser, 

Thr332Ser 

1 

0.57 

±0.04 

(104) 

0.064 

±0.03 

(43) 

0.244 

±0.012 

(19) 

0.01 

±0.008 

(1) 

VI-D2 

Leu69Val,  Hisl  15Val, 
His134Arg, 

Phe222Val, 

Thr332Ser 

2 

1.4 

±0.2 

(255) 

0.33 

±0.2 

(220) 

0.28 

±0.017 

(22) 

0.023 
±0.001 f 
(3) 

MG2-I-A4 

Leu69Val,  Hisl  1 5Ala, 
His134Arg, 

Phe222Met, 

Thr332Ser 

2 

1.95 

±0.4 

(355) 

1.3 

±0.2 

(867) 

1.28 

±0.1 

(98) 

0.13 

±0.006 

(16) 

IV-D11 

Leu69Val,  Hisl  1 5Ala, 
His134Arg, 

Met196Leu, 

Phe222Met, 

Thr332Ser 

2 

2.5 

±0.2 

(455) 

0.73 

±0.3 

(487) 

1.06 

±0.08 

(82) 

0.12 

±0.012 

(16) 
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II-A1 

Leu69Val,  His115Leu, 
Hisl  34Arg, 
Phe222Met, 

Thr332Ser 

2 

2.5 

±0.21 

(455) 

0.73 

±0.5 

(487) 

0.59 

±0.1 

(45) 

0.07 

±0.003f 

(9) 

V-B3 

Leu69Val,  Hisl  1 5Val, 
Hisl  34Arg, 

Phe222Met, 

Thr332Ser 

2 

2.6 

±0.1 

(473) 

0.68 

±0.5 

(453) 

1.2 

±0.5 

(92) 

0.05 

±0.003f 

(6) 

2-II-D1 2 

Lys70Asn, 

Hisl  15Leu, 
His134Arg, 
Phe222Met, 
Thr322Ser 

3 

0.27 

±0.014 

(49) 

0.27 

±0.014 

(180) 

0.96 

±0.15 

(74) 

0.18 

±0.028 

(23) 

1-I-D10 

Leu55lle,  Leu69Val, 
Hisl  15Ala,  His134Arg, 
Phe222Met, 
lle291Phe,  Thr322Ser 

3 

0.87 

±0.076 

(158) 

0.87 

±0.076 

(580) 

2.11 

±0.33 

(162) 

0.31 

±0.063 

(39) 

I-IV-H9 

Leu55lle,  Leu69Val, 
Hisl  15Leu, 
His134Arg, 
Asp136His, 
Phe222Met, 
lle291Leu,  Thr322Ser 

3 

4.2 

±0.5 

(764) 

4.2 

±0.5 

(2800) 

2.38 

±0.18 

(183) 

0.31 

±0.03 

(39) 

VH3 

Leu69Val,  Hisl  1 5Ala, 
His134Arg, 
Phe222Met, 
lle291Leu,  Thr332Ser 

4 

2.6 

±0.2 

(473) 

2.6 

±0.2 

(1733) 

2.8 

±0.2 

(215) 

0.3 

±0.11 

(38) 

VIID2 

Leu55lle,  Leu69Val, 
Hisl  1 5Ala,  His134Arg, 
His197Arg, 
Phe222Met, 
lle291Leu,  Thr332Ser 

4 

3.7 

±1 

(673) 

3.7 

±1 

(2467) 

3.9 

±1 

(300) 

0.17 

±0.014 

(21) 

a.  Mutations  relative  to  the  wild-type  like  rePONI  variant  G3C93.  Mutations  that 
were  newly  introduced  in  a  given  round  are  denoted  in  bold. 


b.  Fold  improvement  relative  to  wt-like  rePONI .  Errors  of  values  were  derived  from 
at  least  two  independent  measurements.  The  maximal  deviation  between 
different  enzyme  preparations  was  <2  fold  (measured  for  VIID2). 

c.  Fast  and  slow  hydrolysis  of  the  two  equally  toxic  isomers  of  GD  (ScSp  and 
RcSp). 

d.  Kinetic  parameters  were  determined  with  the  in  situ  generated  G-agent  and  by 
assaying  residual  AChE  activity.  They  thus  relate  to  the  toxic  Sp  isomer. 

e.  Values  from  ref1 

f.  Standard  deviation 
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Table  7:  Catalytic  activities  of  the  best  variant  from  each  round 


Variant 

Genotype3 

Round 

GDDC 

kcat/Km 

GFD 

kcaf/Km 

GBD 

kCat/Km 

GAD 

^ca/Km 

xIO7 [M 

Fast 

1  min1] 

Slow 

xIO7 

[Mimin'1] 

xIO7 

[Mimin'1] 

xIO7 

[Mimin'1] 

rePONI 

2D8 

Leu69Gly, 
Hisl  1 5T  rp, 
His134Arg, 
Phe222Ser, 
Thr332Ser 

0 

0.0055 

±0.002 

0.4 

±0.1 

0.0015 

±0.0006 

0.015 

±0.003 

0.01 3° 
±0.003 

0.46d 

±0.01 

0.008 

±0.001 

0.025 

±0.002 

0.043 

±0.007 

0.082 

±0.005 

Fold 

improvement 
relative  to 
(re  PON  If 

(73) 

(10) 

(35) 

(3) 

(2) 

PG1 1 

Leu69Val, 

His115Ala, 

His134Arg, 

Phe222Ser, 

Thr332Ser 

1 

1.43 

±0.3 

0.32 

±0.04 

0.2 

±0.03 

0.02 

±0.005 

0.006 

±0.001 

Fold 

improvement  R1 
relative  to 
variant  2D8 
(rePONI  f 

4 

(260) 

21 

(213) 

0.4 

(15) 

0.8 

(3) 

0.07 

(0.1) 

VII-D11 

Leu69Val, 
Hisl  1 5Ala, 
His134Arg, 

Phe222Met, 

Asp309Gly, 

Thr332Ser 

2 

2.9 

±0.9 

2.9 

±0.1 

1.07 

±0.08 

0.12 

±0.01 

0.025 

±0.003 

Fold 

improvement  Ft2 
relative  to 
variant  2D8 
(rePONI  f 

7 

(527) 

193 

(1933) 

2 

(82) 

5 

(15) 

0.3 

(0.6) 

1-I-F11 

Leu55Met, 

Leu69Val, 
Hisl  1 5Ala, 
His134Arg, 

DkAOOO^/lM 

3 

4.4 

±0.6 

4.4 

±0.6 

1.52 

±0.1 

0.39 

±0.04 

0.17 

±0.04 

I  I  V_s  L _ L _ L _  I  V  I  V_s  l  . 

lle291Leu, 

Thr322Ser 
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Fold 

improvement  R3 
relative  to 

11 

293 

3 

16 

2 

variant  2D8 
(rePONI  f 

(800) 

(2933) 

(117) 

(49) 

(4) 

Leu55lle, 
Leu69Val, 
Hisl  1 5Ala, 


IIG1 

Hisl  34Arg, 
Asp136Gln, 

Phe222Met, 
Ile291  Leu, 
Thr332Ser 

4 

5.1 

±0.6 

5.1 

±0.6 

3.4 

±0.3 

0.32 

±0.01 

0.23 

±0.004 

Fold 

improvement  FI 4 

13 

340 

7 

13 

3 

relative  to 
variant  2D8 

(927) 

(3400) 

(262) 

(40) 

(5) 

(rePONI  f 

g.  Mutations  relative  to  the  wild-type  like  rePONI  variant  G3C93.  Newly  introduced 
mutations  are  denoted  in  Bold. 

h.  Kinetic  parameters  relate  to  the  toxic  Sp  isomers  of  GB,  GD  and  GF  or  to  the 
toxic  Rp  isomer  of  GA. 

i.  Fast  and  slow  hydrolysis  of  the  two  equally  toxic  isomers  of  GD  (ScSp  and  Rc 
Sp). 

j.  Values  previously  measured  in1 

k.  Fold  improvement  in  catalytic  activity  of  the  best  variant  of  the  denoted  round 
relative  to  the  starting  variant  2D8,  or  to  wild-type-like  rePONI .  R1  -round  1 ,  R2- 
round  2,  R3-round  3,  R4-round  4. 
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Table  8:  Table  3.  Stereoselectivity  with  Sp-  and  /?p-CMP-coumarin. 


Sp-CMP-coumarin 


Bp-CMP-coumarin 


Round 

# 


Variant 


kcat 

(min'1) 


Km 

(MM) 


kcat/  Km 

(ijM^min'1) 


kcat 

(min'1) 


Km 

(MM) 


l^ca/  Km 
(ijM^min'1) 


Ratio  of 
Sp/  Rp 


1 

2 

3 

4 


rePONI 

<0.0002 

4.5 

54 

(G3C9)a 

n.d. 

n.d. 

(±0.05) 

(±2.6) 

PG1 1 

632 

(±5) 

104 

(±12) 

6.1 

88 

(±2.5) 

381 

(±27) 

VIID11 

500 

105 

4.8 

73 

999 

(±19) 

(±1) 

(±4) 

(±98) 

1-I-F11 

1004 

83 

12.1 

n.d.b 

n.d.b 

(±98) 

(±4) 

VIID2 

1188 

79 

15 

13.8 

1664 

(±24) 

(±2) 

(±1.3) 

(±222) 

~a.  Previously  measured  in' 
b.  n.d.  -  not  detectable. 


0.08 

<0.0025 

0.231 

26 

0.073 

66 

0.003 

4033 

0.008 

1875 

c.  Ratio  of  catalytic  activity  (kcat/KM)  on  Sp-CMP-coumarin  to  the  catalytic  activity  (kcat/KM)  on  RP- 
CMP-coumarin. 
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Table  9:  Progress  in  evolution  of  VX-hydrolyzing  variants 

Catalytic  Activity  (kcat/Km,  M^min1)  of  PON1  Variants  hydrolyzing 
the  toxic  isomers  of  V-type  agents  (inhibition  of  hAChE) 


Variant 

Round 

VX 

RVX 

Wt-G3C9 

<2 

<2 

8C8 

VX-GO 

93 

4E9 

VX-GO 

80 

VIIH3 

VX-R1 

413 

327 

VA4 

VX-R1 

487 

249 

IA12 

VX-R1 

288 

114 

S-2-C8 

VX-R2 

10,600 

<200 

NOTE:  Bacterial  PTE  hydrolyzed  VX  at  8,500  M  imin'1 
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Figures 


Figure  1.  Nerve  agents 

(la)  G-  and  V-type  type  nerve  agents:  tabun  (GA),  sarin  (GB),  soman  (GD),  cyclosarin  (GF)  and 
VX 

(lb)  Coumarin  analogs 
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Figure  2.  Hydrolysis  of  CMP-coumarin  and  CMP-F  by  rePONl  variants 

Enzyme  concentrations  were  varied  depending  on  the  variant's  activity  and  are  noted  in  the 
figure,  (a)  Hydrolysis  of  racemic  CMP-coumarin  (12  pM)  in  the  presence  of  variants  4E9,  3D8, 
3B3  (plus  addition  of  0.03  pM  4E9  after  6  min,  indicated  by  the  black  arrow)  and  wt-like 
rePONl  (plus  addition  of  0.03  pM  4E9  after  20  min);  (b)  Hydrolysis  of  SV-CMP-coumarin  (6 
pM)  in  the  presence  of  variants  4E9,  3D8,  3B3  and  rePONl;  (c)  Residual  AChE  activity  was 
assayed  after  the  incubation  of  in  situ-generated  CMP-F  (40  nM)  and  4E9,  3D8,  3B3  and 
rePON  1  at  the  concentrations  noted  in  the  figure  (data  were  fitted  to  a  first-order  rate  equation  to 
derive  the  apparent  rate  constant  for  hydrolysis  of  CMP-F;  error  bars  represent  the  s.d.  of  two  or 
more  independent  measurements) 
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Figure  3.  G2E6-PON1  penta-mutant  (2D8)  at  2.6  A  resolution  (library  of  round  0).  The  flexible 
loop  (cyan)  covers  the  active  site. 
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Figure  5.  S-2-C8-induced  hydrolysis  of  1.7x10  5  M  VX  using  DTNB  to  capture  the  leaving 
group. 
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Figure  6.  Ex-vivo  stability  of  P0N1  variants  in  human  whole  blood. 

Upper  panel:  assayed  by  monitoring  the  protection  of  blood  ChE  activity  against  GF  inhibition. 
100%  corresponds  to  41-45%  residual  activity  of  blood  ChE. 

Lower  panel:  Evaluated  by  determination  of  PON1  residual  activity  using  the  GF  surrogate 
CMP-coumarin. 
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